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Abstract 

We analyse the Spherical Model with frustration induced by an external gauge 
field. In infinite dimensions, this has been recently mapped onto a problem of q- 
deformed oscillators, whose real parameter q measures the frustration. We find the 
analytic solution of this model by suitably representing the (/-oscillator algebra with 
g-Hermite polynomials. We also present a related Matrix Model which possesses 
the same diagrammatic expansion in the planar approximation. Its interaction 
potential is oscillating at infinity with period log(g), and may lead to interesting 
metastability phenomena beyond the planar approximation. The Spherical Model 
is similarly (/-periodic, but does not exhibit such phenomena: actually its low- 
temperature phase is not glassy and depends smoothly on q. 
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1 Introduction 



In two recent papers |lj, |2|, Parisi et al. have introduced and analysed the Spherical and 
XY Spin Models with frustration, but in the absence of quenched disorder. Their aim was 
to test the conjecture that the frustrated deterministic systems at low temperature behave 
as some suitably chosen spin-glass models with quenched disorder ||. They considered 
the frustrated models in the limit of large dimensionality D of the lattice, where the 
saddle-point approximation becomes exact. In their analysis of these models, they showed 
that the high temperature expansion can be nicely rewritten by using the g-oscillators 
algebra [|T], |2[. Here q measures the frustration per plaquette and varies continuously 
between the fully- frustrated case (q = — 1, fermionic algebra) and the ferromagnetic case 
(q = 1, bosonic algebra). Similar g-deformed algebraic relations have also appeared in 
the Hofstadter problem of quantum particles hopping on a two-dimensional lattice in a 
magnetic field Q. This problem is closely related to the frustrated spin models, which 
can be considered as simplified models of hopping in the large D and classical limits. 
This relation provides another motivation for our analysis. Finally, we note that the 
frustrated XY Model can also describe Josephson junctions arrays in a magnetic field. 

In this paper, we solve exactly the frustrated Spherical Model in the large D limit, 
both in the high and low temperature phases. We use the "coordinate" representation 
of the g-oscillators, which is given by the g-Hermite polynomials || ||. We find that 
the spectrum of the lattice Laplacian is essentially given by a Jacobi theta function; 
thus, it is periodic along the imaginary axis, with period log(g). This g-periodicity does 
not affect the low-temperature phase of the Spherical Model, which is rather standard 
and non-glassy, the effect of frustration being quantitative only. On the other hand, we 
show that the Spherical Model is associated to a Matrix Model |7], which has the 
same diagrammatic expansion in the planar approximation. The potential of this Matrix 
Model is oscillating at infinity, where it has infinitely many minima at approximate 
distance log(g) (see Fig. (2)). Once the corrections to the planar limit are considered, 
its states become metastable by tunnelling to these minima: therefore, we argue that 
this model could have a band spectrum and possibly behave as a spin glass. Similar 
phenomena could also occur in the frustrated XY Model, whose analysis is however left 
to future investigations. 

The models we consider are described by the Hamiltonian 



The complex field <pj 6 C is defined on the sites (labelled by j) of a D-dimensional hyper- 
cubic lattice, and there are nearest-neighbour interactions Uij. Three different models 
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^<j>jUj k <f>k + h.c. . 
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can be obtained by constraining the field as follows: 

N 

(3H G = (3H + J2\<f>j\ 2 , (2) 
(3H S = /3# + ^(l>/ - Nj , (3) 

N 

(3H XY = PH + Y: ^(|0i| 2 -l) • (4) 

The first is the Gaussian Model, which only exists in the high-temperature phase, where 
the mass term dominates the kinetic term ([I]). The second is the Spherical Model ||, 
which uses the Lagrange multiplier \x to enforce the condition J2 |0j| 2 — N, N being 
the total number of sites. This Lagrange multiplier is promoted to a field \ij in the XY 
Model (f|), whose variables are constrained to live on the unit circle \<f>j\ = 1, Vj. 

The couplings Ujk are complex numbers of modulus one and satisfy the relation Ujk = 
Z7L; they are the link variables of an Abelian lattice gauge field, without kinetic term, 
which produces a static external magnetic field. This is chosen in such a way as to give 
the same magnetic flux ±B for any plaquette of the lattice (the product of the four 
U's around the plaquette is e ±lB ). Therefore, the magnetic field should have the same 
projection on all the axes of the lattice, modulo the sign. In order to avoid the choice of 
a preferred direction in the lattice, these signs are chosen randoml)|] Jl|]. 

The ferromagnetic spin interaction is obtained for B = 0, i.e. Ujk = 1. Non-vanishing 
values of B induce a frustration around each plaquette, which is maximal for B = it, 
the fully-frustrated case. Of course all the intermediate, partially frustrated cases, with 
< B < 7i, are also interesting. As explained in Ref. HJ, the B = theory should be 
treated with care, because the large D limit is different in this case, and the normalisation 
factor 1/V2D in (@) should be replaced by 1/(2.D). 

In order to study the previous models, the main difficulty consists in finding the 
spectrum of the lattice Laplacian in the presence of the magnetic field, which is defined 
as: 

2D 

(A/), = £ U jk f k . (5) 
fc=i 

Similarly to the Hofstadter problem [[| , there is a competition between the periodicity due 
to the lattice and those induced by commensurate magnetic fields of the form B = 2irr/s, 
with r, s integers. Thus, we could expect a complex band structure in the spectrum. 
However, there are simplifications due to the large D limit. The authors of Refs. [0, 

* This is a small amount of randomness, of order D(D — l)/2, to be compared to the ~ L D randomness 
present in systems with quenched disorder. 
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have approached the problem via the high-temperature expansion. In the case of the free 
energy of the Gaussian Model (eq. (0)), this is given by a sum over all closed loops, as 
follows: 

oo i / a \ n oo r>2k 

P F = E „™ WW))» = E?rG t . (6) 

n=0, even 71 \V *U J k=Q ^ 

In this expression, the loops C are arranged according to their length n = 2k, and their 
number is M{n). For each loop, the magnetic field yields a weight, which is given by 
the Wilson loop W(C), the path-ordered product of the couplings U along the loop. The 
brackets ( )„ represent the average of this weight over all the Af(n) circuits of length n. 
In the second expression of eq. (|6]), we introduce the notation Gu for the product of the 
multiplicity and the Wilson loop average. 

Each loop encloses a number of plaquettes and receives a weight proportional to 
exp (iBA) , where A is the sum of plaquettes with signs depending on the orientations 
[[[(J. Due to the average over orientations and loops, the quantity G& is a polynomial in 
the variable 

q = cos B , (7) 

of order k(k — l)/2, which is given by the maximal area enclosed by the loop. Parisi [0 
enumerated these diagrams in the large D limit. First, he showed that the same counting 
is given by the Feynman diagrams with 2k external points, which are joined pairwise by 
lines (propagators) intersecting / times, and have assigned the weight q 1 . These diagrams 
also occur in the topological (large N) expansion of Matrix Models @, where the planar 
limit corresponds to no intersections, i.e. to the q = case. Secondly, Parisi found a 
recursion relation for the coefficients of the polynomial Gk{q) - a sort of Wick theorem - 
which can be nicely expressed by the algebra of the g-oscillators a q , a q : 

a q a\ - qa\a q = 1 . (8) 

These operators [II] act on the Hilbert space spanned by the vectors: \m), m = 0, 1 
as follows, 



, . . . , 



where 



a\\m) = \J\m + l} q \m + 1) , 

a q \m) = y/[m] q \m - 1} , a,|0) = , (9) 

1 - a m 

[m] q = -y— y • (10) 

Using the recursion relation, the weighted multiplicities of the diagrams of eq. (|6|) were 
neatly written as an expectation value over the ground state of the g-oscillators Jl|, 0: 

G k (q) = (0\(al + a q r\0) . (11) 
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2 Spectrum of the Laplacian and High- Temperature 
Expansion 

We shall now make expression (|TT| ) more explicit by introducing the coordinate represen- 
tation for the g-oscillators. First note that in the Gaussian Model @ the quantity 
is nothing else then the trace of (2fc)-th power of the frustrated Laplacian, such that we 
can write the general relation 

Tr [/ (A)] = (0|/«)|0) , x q ee a\ + a q . (12) 

The x q coordinate representation, x q \x) = x\x), is given by the so-called continuous 



g-Hermite polynomials || [12|]. These are defined by: 

H n (x) = (x\n) C n , C n = (\n} q ] .y /2 C , (13) 
where the normalization constant Cq is fixed by Hq(x) = 1 and the g-factorial is 

[n) q \ = [n) q [n-l) q ...[l) q , [l] q = [0] q = 1 . (14) 
These polynomials satisfy, of course, a three-term recursion relation in the index n: 

xH n (x) = H n+l (x) + [n] q H n -i(x) , n > 1 . (15) 

Moreover, they obey a g-difference equation J3J[| in their coordinate x, which ranges over 
the interval x G —2/y/l — q, 2/ \J1 — q . A convenient parametrisation is 

2 

x= cos6, 0e[O,7r]. (16) 

yl — q 

More properties of these g-Hermite polynomials can be found in Ref. ||, where they 
are defined as H n (cos9) = (1 — q) n ^ 2 H n (x). The most important property for us is the 
orthogonalizing measure u q (x) [fl, |Bt |L?H: 



/ v q {x) dx H n (x) H m (x) = 5 n>m [n\ q \ , (17) 

J— 2/ y/l — q 



Vq(x) = ^o 5 -v 8ei (^; 



^ l q fl (-l)"^^ 1 )/ 2 S in[(2n + 1)6} 



71 n=0 



sintf (!-? n ) (l-9 n e 2ie ) (l-g"e- 2ie ) , (18) 



^ n=l 



t This g-periodicity will be better discussed later. 
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where &i(z, q) is the first Jacobi theta function. 

Using these results, we can rewrite the general ground-state expectation value in eq. 
(12) as follows: 



Tr [/ (A)] 



dx Vq(x) f(x) . 



(19) 



Therefore, we have shown that the Laplacian is diagonal in the coordinate representation 



of the g-oscillators: it has a continuous spectrum over the interval [—2/^/1 — q, 2/y/\ — q\, 
with eigenvalue density given by v q in eq. fll8f) . 

Let us discuss this spectrum in some interesting limits. For q = 0, we find Uq(x) = 
1/n yi — x 2 /4, which is the Wigner semi-circle law for the Gaussian Hermitean Ma- 
trix Model 0. This results confirms the previous correspondence between the high- 
temperature expansion of the q = Gaussian Model and the planar Feynman diagrams 
without interaction vertices of the Matrix Model. Actually, this correspondence can be 
extended to all values of q (see Section 4). Note also that the q = frustrated model is 
diagrammatically equivalent to the gauge spin glass where the coupling Ujk are random 
quenched variables JIJ. Next, for q = —1, the measure becomes a representation for 
(S(x — 1) + 8(x + l))/2, and we recover the two states of the fermionic algebra ||. The 
limit q — > 1 is singular, owing to the D — > oo peculiarities said before; nevertheless, v\{x) 
becomes the Gaussian distribution, after multiplicative renormalization [Bfl. 

Let us now check that our expression for the Laplacian reproduces the high tempera- 
ture expansion computed in Refs. [[I], ||. The free energy of the Gaussian Model @ can 
be written in terms of the Laplacian as (3F = — \ogZ = Tr[ln(l — (3 A)]. Therefore, the 
internal energy U (/?) is: 



(3U{(3) = R{(3) 



R(z) = Tr 



1- z A 



(20) 



In the last equation, we introduced the resolvent R(z), which will play a major role in the 
subsequent study of the Spherical Model. R(z) is originally well-defined for real z in the 
interval \z\ < z c = (a/1 — q)/2 and is then analytically extended to the whole complex 
plane of z. By using the results of the previous section, we find: 



R(z) 



1 — z x 



where we introduced 



n=0 



K{a) 



-l) n g n ( n+1 )/ 2 



K 



2z 



-i 2n+l 



i - VT^ 



a* 



a 



(21) 



(22) 



The singularities of R(z) in the complex plane are completely determined by those of the 
simpler function K(2z/y/l — q). Actually, the series is very well convergent for \q\ < 1 
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as any Jacobi theta function. Moreover, for | q\ = 1, it is a geometric series, which is still 
convergent because \K\ < 1 for \z\ < z c . 

The high-temperature expansion of the internal energy of the Gaussian Model can 
be obtained by expanding R(f3) in series; there only appear even powers of [3 with q- 
dependent coefficients. We have computed the series to order 0(/3 18 ) with Mathematica 
| fL4}| , and obtained the polynomials Gk(q) , k — 1, . . . , 9 of eq. (||). We have verified that 
they indeed match the results of Refs. p[, which were found by direct enumeration of 
the graphs on a computer. 

For q = 0, we have U(/3) = ((3 — K{2(3)) / (3 2 and the high-temperature expansion 
is singular at (3 C = 1/2, which is the critical temperature of the Gaussian Model. The 
specific heat diverges as ((3 C — (3)~ 1 ^ 2 at the transition, i.e. the value of critical exponent 
a is 1/2. This result still holds for generic — 1 < g < 1 as a consequence of the relation 
between the singularities of R(z) and K(z). 



3 Solution of the Spherical Model 



We are now ready to solve the Spherical Model. Our approach essentially follows the 
original Berlin-Kac solution of the ferromagnetic model in D = 1,2,3 ||. Indeed we 
shall see that the mechanism of the phase transition is the same, only the form of the 
Laplacian is different. The partition function is 



rfiO+ioo r r , 

Z = d/i V<f> exp /3£ 0tA0-/.(V 

J LM—ioo J 



N 



' flQ — lCO 

r^o+ioo 



rfio+ioo 

/ dfi exp [-Tr[mO - f3A)] + fiN] 

J fiQ — ioo 



(23) 



The integration path over /i is along a straight line which runs in the complex plane 
parallel to the imaginary axis; a positive real part has been added in order to make the 
integration over <fi convergent. The trace in the exponent is proportional to N (A is an 
iV-dimensional operator): in the large N limit, we can apply the saddle point method to 
evaluate the integral over /i. The saddle point equation is: 



Tr 



1 



H-/3A 



and must be solved for real, positive values of \i. Introducing the variable z 
equation can be written in terms of the resolvent (|2T)D, as follows: 



z 



R(z) 



(24) 
(3/ Hi this 

(25) 
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In the case of q = 0, this equation can be easily solved: 



4z 2 ) = - 

z 



z{(3) 



p 



(26) 



1 + P 2 ' 

which corresponds to y = 1 + (3 2 . This solution is valid for z < z c = 1/2, namely it 
describes the high-temperature phase (3 < (3 C — 1. Indeed, let us follow the saddle point 
in the z plane, starting from (3 = 0, where z = too. Upon increasing (3, z moves toward 
the critical value z c = 1/2; for (3 > 1, z cannot increase anymore because it finds the 
square-root branch cut of R(z), and therefore "it sticks to the singularity" ||. Hence, 
for (3 > 1, the saddle point equation has no longer an acceptable solution; nevertheless, 
the leading contribution to the /i integral in (|23|) still comes from the neighbourhood of 
z = 1/2. Therefore, the low-temperature solution is z((3) = 1/2 independent of (3, i.e. 
H = 2(3. The free and internal energies are therefore given by 



(3F 
U((3) 



Tr[ln(/i - (3A)} - y, 
(dp + (d^)d,)((3F) 

1 - Tr 



1 



1 



1-zA 



-Tr 



1-zA 




(27) 



(28) 



The last term in this equation vanishes in the high-temperature phase, because the saddle 
point condition is satisfied; however, it should be included for the low temperature regime. 
As a matter of fact, the result can be written in both phases as follows: 

1 1 

More explicitly, 



U({3) 



(29) 



U(T) 



T 



(q = 0) 



(30) 



T > 1 
T-2 T < 1 

The internal energy is therefore continuous at the transition, together with its first deriva- 
tive. But the specific heat presents a discontinuity in its first derivative. These results 
are qualitatively very similar to those of the Berlin-Kac analysis || ; nevertheless, their 
model is quite different, in the sense that frustration is there completely absent. 

The solution for any — l<g<lisa straightforward generalization of the q = case. 
The discussion of the saddle-point equation Q2"5| ) is analogous, because the singularities 
of R(z) are still given by a square-root branch cut. Although we cannot find an explicit 
expression for z((3), we can follow its behaviour: starting from low (3, where z ~ (3, the 
saddle point moves towards higher values of z until it hits the cut of R(z) at z = z c = 
(y/1 — q^j 1 2. In the low-temperature phase, the value of z sticks to z c and thus fi = (3/z c 
for (3 > (3 C . The critical temperature is given by 



(3 C = z c R(z c 



n=0 



■\\ n q n ^ n+1 ^ 2 



(31) 
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The internal energy is still given by the general formula fl29|) , which now reads 
U{T) = { . ^? , „ „ f . x 6 (-K g < 1) . (32) 



A plot of the internal energy is shown in Fig.(P for the three values q = (—1/ V2, 0,1/ y/2). 
This can be easily obtained without solving the saddle-point equation, by using z 6 [0, z c ] 
as a parameter for both £/(/3) in eq. (|29| ) and /3 in eq. (j25|). Note that for q ^ 0, the 
first derivative of the internal energy is not continuous at the transition. The critical 
temperature goes to infinity for q — > 1 and to zero for q — > — 1, respectively. In these 
limiting cases, we can sum the series in R(z) (PH) and obtain the explicit form of z(/3) 
and the internal energy. For q — 1, we find the meaningless results z = f3 and t/(T) = 
in the high-temperature phase, which collapses to a point (T c = oo): this case has a 
pathologic -D — > oo limit, as discussed in Ref. 0, and a sensible theory should include 
1/.D corrections. 

The limit g = — 1 is well defined and corresponds to the fully frustrated model at 
D = oofj]. Here, the high-temperature phase extends down to T = T c = 0, i.e. there is 
no transition. Actually, we find that the resummed R(z) no longer has the square-root 
branch cut of K{z) and that the mapping (3 = (3(z) is invertible on the whole positive (3 
axis: 

P = y^ 2 , (<?=-i) 



U = - z = — \l-^/l + 4p) . (33) 

Let us remark that the absence of a phase transition has also been observed in the fully- 
frustrated long-range XY model |13| . 



In conclusion, the qualitative behaviour of the (infinite-dimensional) frustrated Spher- 
ical Model in the low-temperature phase is rather smooth and standard: indeed, this 
model does not present the potential features of: (i) glassy behaviour (many ground 
states), which is usually found in systems with quenched disorder; (ii) commensurability 
in the spectrum for B = 2irr/s, which was observed in the D = 2 hopping model in 
a magnetic field. Regarding point (i), we would like to remark that the corresponding 
model with quenched disorder, i.e. the Spherical Gauge Glass, does not exhibit a glassy 
phase either. Actually, this model corresponds diagrammatically to the q = frustrated 
model. In intuitive terms, the spin variables are only loosely constrained by the spherical 
condition, J2i \(pi\ 2 — N, and thus can globally adapt themselves to any complex coupling 
configuration. As far as point (ii) is concerned, the commensurability effects are probably 

^Note that the limits D — > oo and q — * — 1 do not commute, as discussed in Ref. Q. 
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washed out by the D — > oo limit, which oversimplifies the geometry of hopping on the 
lattice. 

On the other hand, we expect a glassy phase in the frustrated XY Model 0. In this 
model, the Lagrange multiplier fi becomes a local field fa and the saddle-point equation 
( f24"D is functional. Nothing changes in the high-temperature phase, where the solution 
Hj = [i = const, is correct. However, the solution(s) in the low-temperature phase can 
be rather different from the one of the Spherical Model, and will not be discussed here. 
Rather, we shall approach this problem from a different perspective, by establishing a 
relation with the well developed subject of Matrix Models. 

4 Analogy with the Matrix Models 

Let us first recall the solution of the Hermitean Matrix models in the the large N ap- 
proximation, which corresponds to the planar Feynman diagrams 0. These models 
are characterized by the (N x N) matrix variable M = M' and by the Hamiltonian 

(3H = TrV(M) = Tr (M 2 /2 H ). After diagonalization of the matrix, one is left with 

the partition function over the eigenvalues Aj, i = 1, . . . , N : 




This can be thought of as being the statistical mechanics of N charges with coordinates 
Aj in one dimension, which repel each other logarithmically and are kept together by 
the external potential V(A). This one-dimensional gas of charges usually has a unique 
phase. In the large N limit, the partition sum is dominated by the contribution of 
the saddle point, which corresponds to the equilibrium configuration of the charges, 
neglecting fluctuations. Moreover, the charges become a continuum with density 
which is normalised to one by a convenient rescaling J7|. The saddle point equation is: 

1 V'(\) =P f a dx , A G (-a, a) , /" dx v{x) = 1 , (35) 

2 J— a A — X J— a 

where P stands for the principal value of the integral. This is an equation for the unknown 
v{x), as a function of the given potential V(X). Following Ref. [0], we introduce the 
function 

F(X) = r , (36) 

J -a A — X 

with A taking values in the complex plane. This function is analytic outside the segment 
of the real axis (—a, a) corresponding to the spectrum; moreover, it goes to zero at infinity 
as 1/A, due to the normalisation of v. For A inside the spectrum, we have 

Re F(X) = -V'(A) , Im F(X) = —tiu(X) , A G (—a, a) . (37) 
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The saddle point equation is thus equivalent to these relations for the function F(X): 
they can be usually solved by analyticity arguments, and determine the density u(X). 

Here we would like to remark that these formulas are rather similar to those encoun- 
tered in the Spherical model. Actually, we can identify the two saddle-point equations 
( f25|) and (|35| ) as follows: 

F < A > s i fl (i)- zs \- (38) 

More precisely, the saddle point equation for the Spherical Model is discussed for A outside 
the spectrum, v(X) is given and (3 = (3{z) is the unknown. On the other hand, in the 
Matrix Model A is inside the spectrum, V(X) is given and v(X) is the unknown. 

We can now define the Matrix Model corresponding to the frustrated Spherical Model, 
as the one which possesses the same eigenvalue density u q (X), eq. (0), in the planar limit. 



By using eqs. (|3q) and (p7|) , we can determine its potential: 



V (A) = 2Re F{X) = 2J1 - q ^(-l)" g n(n+1)/2 cosh [(2n + 1) X ] , (39) 



n=0 

where we introduced the convenient parametrisation 



2 

A= n coshx, (40) 

V 1 - Q 



for A outside the cut (— 2/y/l — q, 2/^/1 — q). The Gaussian Matrix Model with V(A) = 
A 2 /2 is indeed recovered for q = 0. 

The relation (38) between the Spherical and Matrix Models can also be understood 
at the level of diagrammatic expansions. In the former model, R(z) generates the high- 
temperature expansion, whose diagrams were cast into the form of n lines joining 2n 
points with weight q for each intersections In the latter model, F{l/z)/z is the 
generating function for the observables (TrM 2n ) = (J2i^f n ) 0- Their diagrams also 
have 2n external points: in the Gaussian Model (q = 0), they are joined by propagator 
lines with no intersections, due to the planar limit. For q ^ 0, the Spherical Model 
diagrams have non-planar intersections, which are reproduced by interaction vertices in 
the planar Matrix Model diagrams; actually, the potential V(A) contains interactions of 
any order. 

The Spherical Model in the high-temperature phase and the Matrix Model in its 
planar limit have corresponding saddle point equations but rather different free energies, 
so they are not quite the same physical problem. The analogy with Matrix Models is, 
nevertheless, interesting because it could provide some useful technology || for solving 
the XY Model. Furthermore, the Matrix Model itself could exhibit interesting physics 
beyond the planar approximation: in the following, we shall put forward some educated 
guesses which are based on the well-known physical picture of the gas of charges. 
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Figure 2: The first derivative of the Matrix Model potential in eq. (39) as a function of 
X, for q = -0.707. 
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The properties of the Matrix Model potential V^(A) can be found by analysing equation 
(BSI). A plot of V'(x) is shown in Fig. <g), for q = -l/y/2: away from the origin, an 



oscillating behaviour sets in with period| x ~" > X ~ (log?) 5 an d amplitude growing to 
infinity This q-periodicity of the potential can be found by inspection of (j39|), and reads: 

V' (fVV) = 2 (^i-^ V2 e* - q^V* V (e*) . (41) 

This is not an exact periodicity due to the presence of the additive term. However, this 
term becomes negligible for x f ar away from the origin, i.e. A > 1, because V grows 
more than exponentially (the amplitude of fluctuations is of order O (exp(2% 2 /| logg|))). 
The g-periodicity of the potential corresponds to a true periodicity of the eigenvalue 
density v q {6) oc Qi(9/n, q), because these two quantities are the real and imaginary parts 
of the function F in (|36|), respectively (their variables (x, 6) in ( |16D and (A, x) in ( ^0| ) 



are also related by analytic continuation). Actually, the eigenvalue density is periodic 
in the direction of the imaginary #-axis: v q {6 — {i/2) logg) ~ v q {0), for < q < 1, and 
v q {6 + tt/2 - (i/2) logg) ~ v q (6), for -1< q < 0. 

This periodicity can be interpreted as the potentiality for metastable states, which, 
however, are not realised in the Spherical Model, owing to its simplified dynamics. Indeed, 
its saddle point equation involves the function F(X) , which is also g-periodic, as its real 
part V; however, F(X) never grows sufficiently high to develop the oscillating behaviour 
and, in fact, goes monotonically to zero at infinity. 

On the other hand, the states might become metastable in the Matrix Model beyond 
the planar approximation. In this approximation, the charges form an equilibrium config- 
uration determined by the minimum of V(X) at the origin, and the tunnelling of charges 
to lower nearby minima is suppressed ||: thus, the g-periodicity of the potential is not 
felt. Beyond this approximation, tunnelling switches on and the states in the spectrum 
can become metastable. Clearly, a detailed analysis is necessary to understand the effect 
of tunnelling into a g-periodic set of local minima separated by ever-rising barriers. This 
might lead to a complex pattern of metastability, which is a characteristic of the glass 
phase. 

In conclusion, there is the possibility that this Matrix Model might describe some of 
the expected effects of frustrated magnetic systems in finite dimension D || || . Moreover, 
it should be solvable by known techniques M. The frustrated XY model is another system 
which might develop these effects in the low-temperature phase; it would be interesting 
to pursue the relation between the XY model and the Matrix Model beyond the planar 
limit. 



The period is | logg|/2 for < q < 1. 
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